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Abstract 

This note is the updated outline of the article "Interpolation^ properties of planar 
spiral curves", Fund, and Applied Math., 7(2001), N2, 441-463, published in Russian. 
The main result establishes boundary regions for spiral and piecewise spiral splines, 
matching given data. The width of such region can serve as the measure of fairness of 
the point set, subjected to interpolation. Application to tolerance control of curvilinear 
profiles is discussed. 

1. Introduction 

This note is intended to present in English the results of the article [1], which has arosed 
interest among researchers, working with spirals, i. e. planar curves with monotone curvature. 

It is not a pure translation of [1]. Some auxiliary results of the article, such as state- 
ments (6), were originally developed for "very short" spiral arcs, namely, one-to-one pro- 
jectable onto the chord. These "interpolational properties" are now essentially expanded, 
and presented to readers as the elements of the theory of spiral arcs (see [3, 4] and references 
herein). Here corresponding proofs are omitted. 

Only final results of article [1] are reproduced in this note. They establish boundary 
regions for spiral and piecewise spiral splines, matching given planar point set. The width 
of such region is considered as the objective measure of fairness of the given data. 

Section 2 and Appendix A, absent in the original article, consider the subject in the view 
of tolerances control of curvilinear profiles. 

2. Demonstration of the bounding region 

We consider problems, related to interpolation of planar point sets. Often one can more 
or less definitely decide: "this interpolant seems to be good, and this one is not" (e.g., 
contains unwanted oscillations). In brief, below we consruct the region, enclosing all "good 
interpolants" . 

To begin with, let us demonstrate such regions in the context of two problems, interpolation 
of function vs curve interpolation. Let points (xj,^) be defined by some function y = f(x) 
as in Fig. 1. For the point set a), among various interpolants, we definitely reject non- 
monotonous ones, because the sequence %ji is monotone. We assume that the designer of this 
point set takes care to present all existing extrema of the function to improve the quality of 
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Fig. 1. Trivial bounding regions for function interpolation 



future interpolation. The region, bounding all monotone interpolants, is shown as the union 
of rectangles. 

Now consider the sequence of non-monotone data with maximum at point 5 (point set in 
Fig. l(b,c)). If the point (x 5 ,y(x 5 )) definitely corresponds to the maximum of y(x), bounding 
region can still be constructed (the case b)). 

For function interpolation these regions are trivial and not interesting. Analogous regions 
for curve interpolation turned out to be much more interesting and useful. Drawing with 
curvilinear profile ABCD in Fig. 2 yields an example. The coordinates (xi, yj) of the inter- 
polation nodes could result from some known curve [x(t),y(t)], which might have an exact 
representation in rather complicated expressions or, e.g., in terms of differential equations. 
In real practice such profile would be presented by 20-50 points. We do with only 4 points, 
just to be able to examine the bounding region without a microscope. 

Data for curve interpolation include values of two functions x(t) and y{t), with no argument 
values tj. Monotonicity and extrema of these functions have no importance: these features 
are not even invariant under rotations. To obtain some invariant conclusions, one would 
like to analyze the behavior of curvature k(t). The analogues of "unwanted extrema" of 
f(x) might be extra extrema of curvature. Contrary to extrema of function, whose minimal 
number is directly visible from the plots like Fig. 1, a minimum of curvature extrema is 
far from being evident. Nevetherless, given data Xi, yi (supplemented with tangents at the 
endpoints) allows us to detect if the original curve could be a spiral. 




Fig. 2. Drawing with curvilinear profile ABCD; blue: bounding region, enclosing all possible spiral 
interpolants; green: bounds for the same profile with 7 interpolation nodes. 
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Fig. 3. Construction of the bounding region for spiral interpolants (simple version) 



Below we prove that all possible spiral interpolants stay within the narrow region (shown 
in blue in the above example ABCD, of the width Conversely, any curve, going 

beyond these limits, is not a spiral. 

3. Geometric preliminaries 

First we explain the simplest construction of the bounding region. Denote 

n((p) = (cos<£>, sin<^) 

the unit vector, associated with the angle ip. A set of 12 points Pi, ... , P i2 and end tangents 
n(ri), n(r 12 ) in Fig. 3(a) is the data to be interpolated. In Fig. 3(b) 12 circular arcs are 
traced as follows: arc PiP 2 of the curvature q% passes through points Pi and P 2 , matching 
at Pi given tangent n(n); circular arcs PiP 2 P 3 , P2P3P4, . . ., P10P11P12, with curvatures 
q 2 , ■ ■ ■ , <7n pass through triples of consecutive points; arc P11P12 of the curvature g i2 matches 
given tangent n(ri 2 ) at the endpoint. Thus we obtain, on each chord PjPj + i, two circular 
arcs, forming the lens. Together these lenses form the sought for region. 
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(4) 



Let n(i?j-) denote the tangent to the circular arc Pj^\PjPj + i at Pj, measured in global 
coordinates. Note that 

Tj =/ij_i+/3j_i = fij + Ctj =>- /3j_x = Pj + Ctj, 

>>J I' J i • 'lj I = /', + V > %-i = Pi + C ; - 

3.2. Spiral arc in normalized position 

We consider spiral and circular arcs between two neighbouring nodes as curves, one-to- 
one projectable onto the subtending chord. Therefore they can be represented as functions 
y(x) = fj(x) in the local coordinate system such that the j-th chord becomes the segment 
[—Cj, Cj] of the local x-axis: 

y" 

y(±Cj) = 0; y (-cj) = tanaj, y (c,) = tan/3,-, Kj(x) = - (5) 

(i + yr) 

is monotone curvature function. 

The circular arc, resting on the segment [— c, c] of the x-axis with the tangent vector n(ip) 
at the start point, is denoted as 
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The spiral segment has tangent n(a) and curvature a at the startpoint A — (— c, 0), and 
tangent n(/3) and curvature 6 at the endpoint 5 = (c, 0). Below two biarcs AJB illustate this 
data. Assuming that the curvature (5) is increasing, a = k(— c) < k(c) = b, the following 
relations are valid: 













^ -r^B 



a + (3 > [or sgn(6— a) = sgn(a+/3)] ; (6a) 
ac < — sin a, sin /3 < be; (6b) 

A (a;; c, — /3) < /(x) < A(x; c, a), \x\ <c. (6c) 

Eq. (6a) is modified Vogt's theorem: (see [3], St. 3). Eqs. (6b) are commented in [3], St. 6. 
Eq. (6c) is lens theorem (see [3], St. 5). These relations turn to equalities if the curvature is 
constant. 




4. Bound for a spiral spline 

Interpolation of a curve usually assumes that the neighbouring nodes are located rather 
closely. The corresponding restrictions could be expressed in terms of sufficiently small 
turning angles pj, e.g., \pj\ <tt/2. We apply even more weak constraints, 

Cj^i+Cj cos pj > and Cj+Cj-i cos pj > 0. (7) 



This requires each of two arcs, Pj-iPj and PjPj + \, of the circle Pj^iPjPj + i not to exceed 
180°. Three last examples below (red) violate these conditions. 




p=33° p=50° p=90° p=95° p=139 d p=160° p=139° p=118° p=96° 

Propositions 1 and 2 constitute Theorem 5 in [1]. 
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Proposition 1. The union of lenses, formed by arcs A(x; Cj, —rjj) and A(x; Cj,^j) on every 

chord PjPj + i, j = l,..., M, covers allpossible spirals, matching given data {Pi, . . . , Pm+i] t i, t ai+i}- 

Proof. In the below proofs only the case of increasing curvature is considered. Denote fj(x) 
the j-th segment of original (or interpolated) spiral. Inequalities near fj(x) in the chain 

A(x; cj, -rjj) < A(x; Cj, -fa) < fj(x) < A(x; Cj, atj) < A(x; Cj,Q (8) 

constitute the lens theorem for a spiral arc (Eq. (6c) and Fig. 4(a)); non-strict form "<" takes 
into account the possibilty of constant curvature. Outer inequalities express the statement 
of Prop. 1 for every j-th chord. To show that the lens, defined by (ay, — j5j), is enclosed by 
the lens —rjj), we have to prove 

> > -ft > -Vr (9) 

The inner inequality, a+(3 > 0, is Vogt's theorem (6a) for increasing (>) or constant (=) 
curvature of the spiral arc fj(x). Note also that, by construction, 

«i=6, 0m=Vm- (10) 

To illustrate proof of (9), the circle of curvature of the spiral at node j (j = 2) is traced 
by dashed line in Fig. 4(c); tangent n(a 2 ) to the spiral at P2 is also tangent to this circle. 
It is known that a spiral arc with increasing curvature intersects the osculating circle from 
right to left. The point Pi is located to the right of the circle, and the point P3 to the left. 
Therefore the circular arc, passing through points Pi,2.3, intersects the osculating circle from 
right to left, and £2 > o>2- The case £2 = «2 occurs if the segment P1P2P3 of the original spiral 
is coincident with the osculating circle. Similarly, from mutual position of the circular arc 
through P2,3,4, and the osculating circle at P 3 (Fig. 4(d)) we deduce r] 2 > (3 2 - □ 



The width of the united region is 



= max 0,-, 7 - = Cj 

l<j<M J J 
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c j \<H-Qj+i\ 



Increasing N decreases Cj and \qj—qj+i\ ~ \kj—kj +1 \, which means cubic convergence of the 
width with N. Fig. 3(c), showing the region for 16 points, illustrates this effect. 

Proposition 2. If turning angles pj obey constraints (7), the sequence qj, j = 1, . . . , N , of 
3-point curvatures is monotone. 

Proof. Comparison of (7) and (3) yields |£| . < |, and \rj\ - < |; so, 

W (3) 

~Vj — ?j =^ —sinrjj < sin£j =^ —qj + \Cj < —qjCj =^ qj + \ > qj. □ 

5. Bound for spline with distinguished vertices (piecewise 
spiral spline) 

Consider the case when the sequence of 3-point curvatures is not monotone, although 
conditions (7) are satisfied. This means that the original curve is not a spiral. Assume that 

nodes with minimal/maximal curvature qj are exactly the vertices of the curve; (Ha) 
there is at least one node between two neighbouring vertices. (lib) 

Under these assumptions there is also a simple way to construct the bounding region. In 
Fig. 5(a) the ellipse is subdivided by 13 points. The 3-point arcs Pj_iPjPj + i are traced only 
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Fig. 5. Constructing bounds for curves with distinguished vertices (dotted line is the original ellipse) 

for the triples, wherein the midpoint P, is not a vertex. As before, we obtain the pair of 
boundary arcs on chords 2-3, 3-4, 6-7, 11-12, 12-13. We have only one circular arc on each 
chord, adjacent to a vertex. Lacking arcs (12) are shown in Fig. 5(b) by bold (red) lines. 
The angles in (12) are such that, e. g., bold arc 1-2 shares tangent at P<i with the arc 2-3-4; 
tangent to the arc 2-3-4 at P4 is matched by the 2-point arc 4-5, and so on. 

The circle of curvature of the ellipse at the point Pg (Fig- 5(a)) shows that the distances 
(P/Psl, IP3P9I are too big, compared to the local radius of curvature. As consequence, the 
estimated bounding regions 7-8 and 8-9 are too wide. 

Proposition 3. Under assumptions (11) the pair of boundary arcs for the point set with 
distinguished vertices is defined as follows: 

vertex at Py. vertex at Pj+i: boundary arcs for the chord PjPj + \: 

yes no A(x;cj,-T}j), A(x; Cj, - Pj)] ,^ 2 \ 

no yes A(x; c j} r) j+1 - p j+l ), A(x;c j ,Q\ 

no no A(x; Cj, —r/j), Aj(x;Cj,(,j). 



Proof. The case no /no is copied from (8). Condition (lib) excludes the case yes/yes. 

First, we note that for the closed curve there is no need to specify the boundary tangents 
Ti, tat: to define the turning angles pi, Pn definitions (1) should be replaced by P = P/v 
and Pjv+i = Pi, assuming M = N (the number of chords is equal to the number of nodes). 

Consider the segment (3-4-5-6-7) with curvature decreasing in (3-4-5) to the vertex at 
node 5 (minimum), and then inreasing in (5-6-7). Inequalities (9) look like 

P4P5: £4 < «4 < -/?4 <? , P 5 Pe: - V5 < -A; < «5 < I , 

and yield —0:5 < r/5, and 04 < —£4. Missed right-most restrictions can now be obtained as 
follows: 

(4) 

"At = -"5 - P5 < ??5 - P5] "5 = Pa ~ P5 < ~U ~ P5- 

Similar reasoning yields similar bounds for the vertex of maximal curvature; e. g., for vertex 8 
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Fig. 6. Simple bounding region (clashed), and narrowed bounding region (solid lines); biarc boundaries 

are shown with an arrow at the joint. 

Together they look like (with upper signs for a vertex-minimum) 

s, $ <\j -Pj § Vj+i - pj+i ; -Vj £ -Pj $ oij § i j i - pj . 

V v ' V v ' 

chord PjPj+i, vertex at Pj+i chord PjP,+i, vertex at Pj 

This yields bounds (12) for yes/no and no/yes cases. □ 

Note that the discrete curvature plot may look like including the 

segments of constant curvature, which can be named "extended vertices". Any node of the 
extended vertex can be chosen as the vertex to apply the above described algorithm. 



6. Narrowing the bounding region for spiral splines 

Fig. 6 illustrates the strengthened version of Prop. 1: some of the bounding arcs are re- 
placed by biarcs. As the consequence, the bound for the whole region turned into a pair of 
smooth curves, intersecting at the nodes. 

If 4 points Pj-i, ■ ■ ■ , Pj+2 are cocircular (e.g., collinear), we obtain qj = qj+i, and zero 
width 0j = O. Any spiral, matching such data, includes arc Pj^iPjPj + iPj +2 of constant 
curvature qj. "Simple" construction, defined by (8), yields the region of zero width on 
the segment PjPj + \. Strengthened version (18) returns zero width on the whole segment 
Pj_iPjPj + iPj +2 Example in Fig. 7 illustrates this situation. 

To describe the narrowed region (Prop. 5) we first introduce notation for biarc curves. 
Because all involved biarcs are one-to-one projectable onto the chord, they can be considered 
as functions B(x; . . .) in the local coordinate system (6). The first circular arc A J of a biarc 
in (6) has tangent n(a) and curvature a at the startpoint A = (— c, 0), and is smoothly 
continued at the join point J by the second arc JB of curvature b to the endpoint (c, 0) with 
end tangent n(/3). The condition of tangency of arcs A J and JB looks like 

(ac + sin a) (be — sin (3) + sin 2 uj = 0, where to = (13) 

is the angular half-width of the associated lens. Details and reference formulae for a family 
of biarcs are given in [4] . 
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Simple bounding region for the above data 
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Fig. 7. Regions for data with 4 cocircular (collinear) points 




Fig. 8. Defining "degenerate" biarcs (15a) and (15b); lens boundaries are shown dashed 

It is well known that, fitting two-point G 1 data (tangents at the endpoints) with a biarc, 
we have one degree of freedom: one curvature, either a or b, can be selected, the other should 
be defined from condition of tangency (13). We use the following notation for biarc curves: 

Bi(x; c, a, (3, a) for biarc, matching end tangents n(a), n(/3), and the start curvature a; 
B 2 (x; c, a, (3, b) for a biarc, matching end tangents n(a), n(/3), and the end curvature b; 
B (x; c, a, (3,p), p G [0; oo], for a family of short 1 biarcs with curvatures 

be — sin (3 



- (sin a sinw], b = - (sin f3 + psinu) 

c \ V J c 



P 



ac + sin a 



sinu 



Now define degenerate biarcs as follows. For a + (3 ^ 0, and either p = or p = oo: 
B (x; c, a, (3, 0) = lim B Q (x; c, a, (3,p) = BAx; c, a, 0, =Foo) = A(x; c, —0)\ 

B (x; c, a, (3, oo) = lim B (x; c, a, (3,p) = B 2 (x; c, a, (3, ±oo) = A(x; c, a). 



(14) 



(15a) 



1 The whole family includes also long biarcs with p < (see [4]). 
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For a + = 0: 



B Q (x; c,7, -7,p) = limS (x; c,w+7,o;-7,p) = A(ar; 0,7). 



(15b) 



a p 

Fig. 8(a) illustrates definitions (15a): when, e.g., the join point tends to startpoint A, the 
first curvature a tends to infinite impulse of curvature at A, the first subarc of a biarc 
vanishes; simultaneously the second arc tends to the lens boundary. Similar definitions for 
decreasing curvature can be written with a+f3 < (6a), and either a — > +00, or b — > —00. 

In Fig. 8(b) the case a+(3 = 2u — > (15b) is illustrated with lenses, narrowing down in 
width. 



Proposition 4. Consider a spiral arc y = f(x) with increasing curvature n(x). If boundary 
angles a, (3 and curvature fall in ranges 

a' <a< a", P' < f3 < 0", -00 < a < k(x) <b< +00, 
than the curve is bounded by 

ifa'+P" > 0: £i(x;c,a',/3",a) <f(x); 
ifa"+/3' > 0: f(x) < B 2 {x;c,a",p',b). 

The proof ([1], Th. 4) includes inequalities (labels correspond to Fig. 9): 

AJ 3 B 



(16) 
(17) 



A.hB 



AJ 2 B 



Bi(x; c, a', (3", a) < Bi(x; c, a, (3", a) < Bi(x; c, a, (3, a) < f(x); 
f(x) < B 2 (x; c, a, (3, b) < B 2 (x; c, a", (3, b) < B 2 (x; c, a", (3', b) 

V v ' V v ' V v ' 

AJ4.B AJ 5 B AJ 6 B 

Inequalities near f(x) result from bilens theorem ([1, Th. 3], [4, Th. 1]). The rest can be 
obtained by inspecting the involved circular arcs. 

Proposition 5. The narrowed bounding region for spiral data {(xj, yj),j = 1, . . . , M+l, Tx, ta/+i} 
with increasing curvature is defined by 

B (x; c v /?>;.) < f 3 {x) < B (x; Cj , afj, /^pj)- (18a) 
Biarcs 7 parameters are 



a" 

j 



aj 



Pi 



3 = 1 



6 



-00 



1 < j < M 



max(-p i -^_ 1 , -rjj) 







c 7 -i 



sin f3';_ x 



sinw 



ajCj + sin a 



1 ■ 



to 



4 + % 



0'j 



P 



l<j<M 



max(-^-, p j+ i-T] j+1 ) 



sm a 



C3+1 



3+1 



j = M 



Vm 



Vm 



+00 



bjCj — sin (3'j 



00 — 



smw 



(18b) 
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Proof. Transform (9), following the scheme: 

-£i < 0i < Vi , -£i~P2 < 0i~P2 < rji-p 2 -£i~P2 <a 2 < £2, 

-7/2 < a 2 < 6 -V2+P2 < a 2 +p 2 < -V2+P2 <Pi<r)i- 

Combining the first and the last column yields ranges (16) for 0i and a 2 . 

max(-^ i ,p 2 -^2) y <Pi< m , max(-p 2 -£n -^) y < «2 < • 

' J[ ' 0'[ ' ^ ' a 2 ' 

Similarly, we obtain ranges for pairs (0 2 , 013), (03, a^), . . . , (0m-i,<xm)- Missing ranges for 
q;i and /3m are known from (10): a[ = a" = £1, /3 M = /3 M = t\m- Two rows of the table (18b) 
are thus filled. 

Now check conditions in (17). E.g., for a^+/3'' > we have 

(9) 

j = l: a ; +/ ^ = £ 1 + 7h > 0; 

j > 1: aJ.+^J-' = max(-7/j, -pj-^^) + rjj = max(0, Vj-pj~^j-i) > 0. 

Similarly, conditions ct"+/3j > hold. 

To estimate ranges for unknown curvatures fcj in nodes j = 1, . . . , M+l, we make use of 
inequalities (6b); non-strict form accounts for the possible case of constant curvature: 

on the chord P\P 2 - k\C\ < — sinai, k 2 c\ > sm0i, 

P 2 P 3 : k 2 c 2 < -sina 2 , k 3 c 2 > sm0 2 ; 



(19) 



Pm-\Pm- k M -iC M -i < -sina M -i, k M c M -i > sm(3 M -u 
last chord P m Pm+i- k M CM<—sina M , k M +iCM > sin/3 M . 



This results in 

sin sin , -sin a,- -sina' 
aj = — < — < kj < < = bj, 

Cj-i Cj-i Cj Cj 

S v ' S V ' 

except 3 = 1 except last node 3 = M+l 

i.e. — 00 < Kx(x) < b\, a 2 < k 2 (x) <b 2 , clm < acjw-(x) < + oo. Prop. 4 can now be 

applied to define boundary biarcs: 

Si(ar; Cj, a'j, 0", a 5 ) < < £ 2 (x; 9, a", /^-, fy). (20) 

It is easy to verify that degenerate cases, arising in (20) with Oi = —00, or b M = 00, or 
a+(3 = 0, result in "simple" bounds (8). E. g., if ot'j+0" = is the case, a'j = —0" = —r]j, and 

B x (x; Cj, a'j, 0"dj) = B (x; Cj, -rjj, r)j,p) ( = b) A(x; Cj, -rjj). 
Eq. (20) is rewritten in (18) in terms of "universal" function Bq(x; a, 0,p) (14) 2 . □ 

Note also that possible continuation of (19) is 

sin£,_i sin/3,_i , sin a,- sinn,- 
qj-i = — < — < kj < < = q j+1 : 

Cj_i Cj_i Cj Cj 

unknown curvatures kj are limited by known 3-point neighbouring curvatures qj±i- The above 
chosen ranges [dj] bj] yield more narrow region. 

It may also happen that limits ±00 for end curvatures (18b) can be specified from some 
additional considerations. E.g., in the 3-point data (M = 2), shown below, one can assume 



2 In (20) "hidden" degenerate biarcs may occur. E.g., the biarc B^{x; c, a", ft, b), regular at first sight, 
becomes degenerate if be — sin/3' =0, and equal to B\{x; c, a", ft , 00) = Bq(x; c, a", ft, 0). This happens 
in the data like shown in Fig. 7, chord 2. 



11 



positivity of curvature, and replace a\ = — oo by a\ = 0. This replaces the lower boundary 
arc A(x; Cj, —rji) (dashed) by the regular biarc Bi(x; Ci, £1, r]i, 0): 




* * * 

The specific occasion of the research [1] was eventuated in about 1995 3 , in the workshop 
of the Institute for High Energy Physics (IHEP, Protvino). Research originated from a 
practical problem of inspection of a certain cam profile: there was, on a small part of the 
360°-profile, a big discrepancy between nominal and measured curves. The reason, expressed 
in present-day terms, was that the design of the profile included very unfair point subsets. 

While exploring the situation, the author has become interested to find some quantitative 
measure of such drawbacks in drawings, possibly, a tool for an expert-metrologist. Purely 
geometrically the problem was shaped as follows: to describe a set of curves with as few 
vertices as possible, matching given interpolation data. The simplest case (spiral, zero ver- 
tices) gave a well restricted solution, justifying such problem setting. Proposition 3 is the 
step towards non-spiral curves. 

A. Application to computer-aided tolerancing 

A.l. Bounding regions in the view of tolerance control 

Bounding region is not quite a new concept for industrial design. Returning to Fig. 2, 
consider the hole, dimensioned as 26±0.2. This can be treated as the bounding region (in 
the form of ring) for the circular profile. The given tolerance admits inaccuracy of order 0.2 
in coordinates of the center of the hole. If declared as the datum element, the hole defines the 
origin of the coordinate system; consequently, any other dimension of the type "coordinate" 
on this drawing should be agreed with this inaccuracy. E.g., setting tolerance of order 0.1 
for the dimension 100 would be an evident error in the drawing. 

Normally such error will be corrected by an expert-metrologist before starting the man- 
ufacturing process. Possible correction could be specifying the roundness tolerance e <C 0.1 
for the hole. This allows the value of diameter still to be in the range 25.8 < < 26.2, but 
reduces the width of the ring-shaped bounding region to much smaller acceptable value. 

Note that, given all dimensions and tolerances, one can usually construct the exact model 
of the ideal part, corresponging, for the case of the hole, to 26±0. The position of any 
point of the circular and straight line elements will be perfectly determined. But this is not 
so for the elements like profile ABCD. Since it is not presented by exact equations, and 
the interpolation method is not exactly specified, there occurs certain nondeterminancy is 
constructing the nominal profile itself. 

Assuming a coordinate inspection machine as a device for tolerance control, one should 
assure precision of coordinate measurements about 10% of tolerances under inspection. The 

3 Article [1] was first published in 1998 as Preprint IHEP, 98-9, Protvino, 1998. 
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Fig. 10. Interpolating data with a cubic spline 



more accurate is the device, the more precisely the distance from the measured point to a 
circle/line can be defined. Again, this is not the case for the curvilinear profile: one could 
define exact distance, say, to the currently chosen interpolant. The latter could be different 
from one, previously used in NC-machining process. 

So, the nondeterminancy of the profile behaves as inevitable addition to the measurement 
error. Constructing the bounding region estimates this nondeterminancy, and allows us to 
set its acceptable limits, called by the accuracy requirements. The width of the region is 
the objective measure of fairness of the given point set, and may serve as the criterion to 
demand a designer for more detailed profile description. 

For traditional interpolation methods nothing ensures that some interpolant passes within 
the prescribed bounding region. But, on the other hand, one can verify whether it does or 
not. An example is discussed below. 
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A. 2. Test with cubic spline interpolation 

Algorithms of interpolation with spirals are not widely distributed with CAD software. 
On the other hand, interpolating spiral data by, e.g., a cubic spline, will usually not return 
a spiral. However, as soon as the region is known, wherein the sought for curve definitely 
lies, acceptance test for any interpolation method can be performed. We add to this note 
an example of such test. In Fig. 10 given data points and end tangents were interpolated by 
a cubic spline curve. The nodes in the example are spaced out such as one could visually 
estimate if the spline fits the prescribed region (three last chords are shown magnified). 

The spline can be parametrized by accumulated chord length, which approximates arc 
length of the interpolated curve. Assuming that the 3-rd order Bezier polynomials x(t), y(t) 
are close to the Taylor series x(s), y(s) of the naturally parametrized unknown curve gives 
rise to set boundary conditions at the start point of the cubic spline as (x r (0),y'(0)) ~ n(rj), 
and similarly at the endpoint. 

Points setting of 5 first nodes, compared to expected curvatures, is evidently unacceptable 
for good interpolation. It is exhibited in both location of the spline out of prescribed region, 
and large beatings in the curvature plot. 

Situation on the right side looks much better, although points are set still too rare, com- 
pared with typical industrial drawings. Recall that the bounding region encloses some 
unknown, but well defined curve (e.g., solution of a differential equation, or some other 
implementation of designer's concept). 

• First, we estimate whether the bounding region is sufficiently narrow, and the nonde- 
terminancy of the profile is possibly negligible. 

• Second, we estimate closeness of the (unknown) curve and the interpolant, which may 
even fall inside the bounds. 

Appearance of the curvature plot in Fig. 10, with curvature extrema in every node and 
between them, is rather similar to the examples in [2, Figures 9.8-9.13]. The question arises 
if the non-monotonicity of the curvature plot is essential, and should be somehow corrected. 

The general answer depends on the values of curvature, and the particular reason, by 
which the beatings are desired to be avoided. Analyzing the problem, one should not for- 
get that the objects of computer-aided design become the objects of further manufacturing, 
which is not geometrically precise. We can compare the model and the real part (cam profile, 
highway path, etc.), and find that the latter fits well prescribed tolerances. But this small 
difference could result in rather big distortions of curvature plot: 



C 




In the above illustration the data is chosen on the circular arc of radius 10 with P\ = (0, 0), 
and 7~i = 0. The coordinates were rounded to two decimal digits: P 2 ~ (0.5233, 0.0137) was 
repalced by (0.52, 0.01), P 3 « (1.0453, 0.0548) by (1.05, 0.05), . . . , and P 21 « (8.6603, 5.0) 
by (8.66, 5.0). This rounding caused well visible beatings in discrete curvature plot. The 
geometric reason of such behavior is evident. For a small circular arc ABC a slight shift of 
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the midpoint B effects a dramatic change in radius or curvature. The smaller is the angllar 
measure of the arc, the greater is this effect. 

It seems that the above considerations agree with the note in [2, Sec. 23.1]: the definition 
of fairness, based on minimum of curvature extrama, "is certainly subjective; however, it has 
proved to be a practical concept". 
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